Time-reversal violating Schiff moment of ^^^Ra 
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We use the Skyrme-Hartree-Fock method, allowing all symmetries to be broken, to calculate the 
time- reversal-violating nuclear SchiflF moment (which induces atomic electric dipole moments) in the 
octupole-deformed nucleus ^^^Ra. Our calculation includes several eflFects neglected in earlier work, 
including self consistency and polarization of the core by the last nucleon. We confirm that the 
SchiflF moment is large compared to those of reflection-symmetric nuclei, though ours is generally a 
few times smaller than recent estimates. 



I. INTRODUCTION 



Experiments with K and B mesons indicate that time- 
reversal invariance (T) is violated through phases in the 
Cabibbo-Kobayashi-Maskawa matrix that affect weak in- 
teractions . The suspicion that extra-Standard-Model 
physics, e.g. supersymmetry, also violates T has moti- 
vated a different kind of experiment: measuring the elec- 
tric dipole moments (EDMs) of the neutron and of atoms. 
Because any such dipole moment must be proportional 
to the expectation value of the T-odd spin operator, it 
can only exist when T (and parity) is violated 0, |1] ■ So 
far the experiments have seen no dipole moments, but 
they continue to improve and even null results are use- 
ful, since they seriously constrain new physics. Whatever 
the experimental situation in the future, therefore, it is 
important to determine theoretically what the presence 
or absence of EDMs at a given level implies about T- 
violating interactions at elementary-particle scales. Our 
focus here is atoms, which for some sources of T viola- 
tion currently provide limits as good or better than the 
neutron Q. 

One way an atom can develop an EDM is through 
T and P violation in its nucleus. Let us assume that 
given a fundamental source of the broken symmetry one 
can use effective-field theory and QCD to calculate the 
strength of the resulting T-violating nucleon-pion inter- 
action. One then needs to connect the strength of that 
interaction to the resulting nuclear "Schiff moment" , 
which, because the nuclear EDM is screened 0, is the 
quantity responsible for inducing an EDM in electrons 
orbiting the nucleus. The Schiff moment is defined clas- 



sically as a kind of radially weighted dipole moment: 

Jd'rpMr){r'-l'^?jr, (1) 

where pch is the nuclear charge density and r^j^ is the 
mean-square charge radius. Recent papers 0,0 have ar- 
gued that because of their asymmetric shapes, octupole- 
deformed nuclei in the light-actinide region should have 
collective Schiff moments that are 100 to 1000 times 
larger than the Schiff moment in -'^^^Hg, the system with 
the best experimental limit on its atomic EDM (J] . Ref. 
suggested that certain many-body effects may make 
the enhancement a bit less than that. The degree of en- 
hancement is important because several experiments in 
the light actinides are contemplated, planned, or under- 
way [gl llOjl. They may see nonzero EDMs, and even if 
they don't we will need to be able to compare their limits 
on fundamental sources of T violation to those of exper- 
iments in other isotopes. 

Perhaps the most attractive octupole-deformed nu- 
cleus for an experiment is ^^^Ra. Though radioactive, 
it has a ground-state angular momentum J = 1/2, which 
minimizes the effect of stray quadrupole electric fields in 
an experiment to measure a dipole moment^. In addi- 
tion, the associated atom has close-lying electronic levels 
of opposite parity and is relatively easy to trap and ma- 
nipulate. As a result, at least one group is at work on a 
measurement in ^^^Ra Here we calculate its Schiff 
moment, attempting to incorporate the effects discussed 
in Ref. j^] through a symmetry-unrestricted mean-field 
calculation. We begin in the next section by describing 
the physics of the Schiff moment in octupole-deformed 
nuclei, briefly reviewing prior work in the process. In 
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^ The statement that the nucleus has octupole and quadrupole de- 
formation really refers to its intrinsic state, a concept we elabo- 
rate on below, and does not contradict its insensitivity to applied 
electric fields with multipolarity greater than one. 
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Section IIIII we test our mean- field approach by calculat- 
ing properties of even Ra isotopes. In Section IIVI we 
discuss issues peculiar to mean-field calculations in odd 
nuclei and then present our results for the SchifF moment 
of ^^^Ra, focusing particularly on the degree of enhance- 
ment. Sectionals a brief conclusion. 



where |^'o) is the member of the ground-state multiplet 
with = J 7^ 0, the sum is over excited states, and Sz 
is the operator 



II. ENHANCEMENT OF SCHIFF MOMENTS IN 
OCTUPOLE-DEFORMED NUCLEI - PREVIOUS 
WORK 

In analogy with dipole moments in atoms, static SchifF 
moments in nuclei can exist only if T is broken. Because 
T-violating forces are much weaker than the strong in- 
teraction, the Schiff moment can be accurately evaluated 
through first-order perturbation theory as 



Eo — Ei 



■ c.c. 



(2) 



— — 

— in 



10 5Z 

P 



^/i-.2 

> 3' 



ch 



(3) 



with the sum here over protons. The operator VpT is 
the T- (and parity-) violating nucleon-nucleon interac- 
tion mediated by the pion [ill ll^ (shown to be more 
important than other mesons in Ref. 
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where arrows denote isovector operators, Tz is -|-1 for neu- 
trons, mjv is the nucleon mass, and we are using the con- 
vention h = c = 1. The g's are the unknown isoscalar, 
isovector, and isotensor T-violating pion- nucleon cou- 
plings, and g is the usual strong TriVA^ coupling. 

In a nucleus such as ^^^Hg, with no intrinsic octupole 
deformation, many intermediate states contribute to the 
sum in Eq. |2Jl. By contrast, the asymmetric shape of 
^^'"^Ra implies the existence of a very low-energy |l/2^) 
state, in this case 55 keV above the ground state |5'o) = 
|l/2+), that dominates the sum because of the corre- 
sponding small denominator. To very good approxima- 
tion, then, 

(l/2+|g,|l/2-)(l/2-|VpT|l/2+) , _ 

^ +C.C., (5) 

where /S.E = 55 keV. The small denominator is part of the 
reason for the enhancement of the Schiff moment. The 
other part is the matrix element of the Schiff operator in 
Eq. In the limit that the deformation is rigid, the 
ground state and first excited state in octupole-deformed 
nuclei are partners in a parity doublet, i.e., projections 
onto good parity and angular momentum of the same 
"intrinsic state" that represents the wave function of the 
nucleus in its own body-fixed frame. The matrix elements 



in Eq. (jSJ are then proportional (again, in the limit of 
rigid deformation) to intrinsic-state expectation values, 
so that [6j 
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where J is the ground-state angular momentum, equal 
to 1/2 for ^^^Ra, and the brackets indicate expectation 
values in the intrinsic state. The intrinsic-state expecta- 
tion value (5'^) is generated by the collective quadrupole 
and octupole deformation of the entire nucleus; it is 
much larger than a typical matrix element in a spher- 
ical or symmetrically deformed nucleus. Together with 
the small energy denominator, this large matrix element 
is responsible for the enhancement of laboratory-frame 
Schiff moments in nuclei such as ^^^Ra. 

The amount of the enhancement is not easy to calcu- 
late accurately, however. The reason is that the matrix 
element of the two-body spin-dependent operator VpT 
in Eq. ^ depends sensitively on the behavior of a few 
valence particles, which carry most of the spin. In the 
approximation that particles (or quasiparticles) move in 
independent orbits generated by a mean field, the po- 
tential can be written as an effective density-dependent 
one-body operator that we will denote Upx, defined im- 
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plicitly by 

{a\l}pT\h)^Y.^ac\VpT\hc), (7) 

c<F 

where |a), |6), and |c) are eigenstates of the mean field 
and the matrix elements of VpT are antisymmetrized. 
With the further approximation that the mass of the pion 
is very large, UpT can be written as a local operator, in a 
form we display in the Section Hvl Evaluating its matrix 
element is tricky. 

The authors of Refs. 0,01 used a version of the particle- 
rotor model Q to represent the odd-A nucleus. In this 
model, all but one of the nucleons are treated as a rigid 
core, and the last valence nucleon occupies a deformed 
single-particle orbit, obtained by solving a Schrodinger 
equation for a Nilsson or deformed Wood-Saxon poten- 
tial. The model implies that the core carries no intrinsic 
spin whatever, that the neutron and proton densities are 
proportional, and that the exchange terms on the right- 
hand side of Eq. Q are negligible. Under these assump- 
tions, UpT, which now acts only on the single valence 
nucleon, reduces to 

UpT{r)^r^ ^^^^^ rT-Vpo{r). (8) 

where G is the Fermi constant, inserted to follow con- 
vention, and po is the total nuclear mass density. The 
dimensionless parameter rj is then a function of the cou- 
plings Oi and the isospin of the nucleus. 

Ref. jPl confirmed the collectivity of the intrinsic Schiff 
moments obtained in Refs. ^Jl, but questioned the ac- 
curacy of some of the assumptions used to evaluate the 
matrix element of Vpr, suggesting that either core-spin 
polarization or self-consistency in the nuclear wave func- 
tion might reduce laboratory Schiff moments. The zero- 
range approximation and the neglect of exchange in UpT 
are also open to question. As a result, it is not clear 
whether the Schiff moment of ^^^Ra is 1000 times that of 
is^Hg or 100 times, or even less. In what follows, we pro- 
vide a (tentative) answer by moving beyond the particle- 
rotor model. Our calculation is not the final word on 
Schiff moments in octupolc-dcformcd nuclei — we only 
do mean-field theory, neglecting in particular to project 
onto states with good parity, and do not fully account for 
the pion's nonzero range — but is a major step forward. 



III. MEAN-FIELD CALCULATIONS FOR 
OTHER Ra ISOTOPES 

A. Mean- field calculations 

Self-consistent mean-field theory is widely used for 
describing bulk properties of nuclei [l6j . In the guise 
of density-functional theory, it is also used throughout 
atomic and molecular physics. The approach is more 



"microscopic" — nucleons are the only degrees of free- 
dom — and far less phenomenological than the col- 
lective particle-rotor model. Self-consistency connects 
the single-particle states and the actual density distri- 
bution. The variational principle that determines the 
single-particle wave functions thus optimizes all multi- 
pole moments not fixed by global symmetries. The den- 
sity distributions of neutrons and protons are not propor- 
tional to each other; they have slightly different deforma- 
tions and radial profiles. In odd-A nuclei, self-consistent 
calculations include rearrangement due to the unpaired 
particle. Rearrangement causes polarization of the even- 
A core through orbital-current and spin-density terms in 
the effective interaction. Core polarization is one of the 
effects on the Schiff moment of ^^^Ra that we investigate 
below. 

Our approach is nonrelativistic and employs Skyrme 
interactions. To get an idea of the range of results this 
kind of calculation can produce, we use four different 
parameterizations of the Skyrme energy functional, i.e., 
four different Skyrme forces. The four give similar re- 
sults for many observables near stability, but still have 
significant differences. Our favorite interaction, for rea- 
sons explained below, is SkO' 0,0], but we also show 
results for the commonly used forces SIII ^ , SkM* , 
and SLy4 



B. Related Observables in Even Isotopes 

Intrinsic-parity breaking in even radium isotopes is the 
subject of several theoretical analyses; see the review in 
Ref. m and the more recent studies in Refs. [2ll2^. To 
assess the ability of the Skyrme interactions to handle it, 
we perform a series of Hartree-Fock (HE) -I- BCS calcu- 
lations for even Radium isotopes. We use the Skyrme- 
HE+BCS code from Ref. 25]; it represents single-particle 
wave functions on an axially symmetric mesh, and uses 
Eourier definitions of the derivative, 1/r^, and l/r^^ op- 
erators. We choose 75 grid points in the z direction, and 
27 in the rj^ (perpendicular) direction, with 0.8 fm be- 
tween them. The code uses a density-independent zero- 
range pairing interaction with a self-adjusting cutoff as 
described in Ref. l26|. Eor each Skyrme force we adjust 
the pairing strength separately for protons and neutrons 
'26^ . We should note that other self-consistent mean-field 
models, namely HE-I-BCS with the nonrelativistic Gogny 
force [23 and the relativistic mean-field model , yield 
results that are similar to those we describe now. 

Figure n illustrates the calculated evolution of intrin- 
sic deformation with increasing neutron number in the 
Radium isotopes. It plots the intrinsic ground-state 
mass-density contours predicted by SkO'. The mean- 
field ground states go from having a spherical shape 
at the magic number N = 126 to a quadrupole de- 
formed (refiection-symmetric) shape at TV = 130, then to 
quadrupole-foctupole deformed (refiection-asymmetric) 
shapes for 132 < N < 140, and finally back to quadrupole 




deformed shapes at higher N. Because the ground 
states are obtained from a variational principle, all shape 
moments higher than octupole are also optimized (the 
isoscalar dipole moment is constrained to be zero). The 
nucleus ^^^Ra, with N = 137, will clearly be well de- 
formed in both the quadrupole and octupole coordinates. 
The structures at small radii visible for N > 132 reflect 
small oscillations of the density distribution around the 
saturation value (for a given neutron excess) caused by 
shell effects. 




124 128 132 136 140 144 148 
Neutron Number N 

FIG. 2: Relative error in binding energy (top) and predicted 
two-neutron separation energies (bottom) for four Skyrme in- 
teractions in a series of even- A'' Radium isotopes. The exper- 
imental separation energies are also shown. 

We must note that the octupole-deformed minima are 



not equally pronounced for all forces and isotopes. In ad- 
dition, in some of the isotopes with reflection-symmetric 
minima, some of the Skyrme forces predict an excited 
octupole-deformed minimum separated by a small bar- 
rier from the ground-state minimum. Furthermore, in 
the transitional nuclei, which have soft potentials in the 
octupole direction, all parity-breaking intrinsic deforma- 
tions are subject to collective correlations as discussed 
in Ref. 0|. The influence of correlations will be small- 
est for the nuclides with the most pronounced octupole- 
deformed minima, usually ^^^Ra and ^^^Ra. This fact 
supports our belief that our mean-field calculations sup- 
ply a good approximation to the intrinsic structure in 
225Ra. 

Figure 121 shows the relative error in the predicted bind- 
ing energies SE = (i^caic--Eexpt)/£^cxpt for all four forces, 
and the predicted two-neutron separation energies, along 
with the measured values. All the forces do a good job 
with binding, which is not surprising given the way their 
parameters were fit. The fact that the error in binding 
for SkO' is nearly constant with iV for > 130 is re- 
flected in the near perfect agreement in the bottom panel 
with the measured two-neutron separation energies S2n- 
The errors in predicted values of S2n around N = 128 
probably reflect the deficiencies of mean-field models in 
transitional nuclei. 

Figure 13 shows three parity- violating intrinsic quanti- 
ties. In the top panel is the ground-state octupole defor- 
mation (33 ^ 4n(r^Y3o)/{3AR^) (where R = 1.2A^/^), as 
a function of neutron number. The trend mirrors that 
in the density profiles shown earlier. At = 136, one 
less neutron than in ^^^Ra, the forces all predict almost 
identical octupole deformation, a result we like. Experi- 
mental data for octupole moments are still sparse in this 
region; we are only aware of = 0.105(4) for N — 138, 
a value that can be deduced from the B{E3;Q+ 3j") 
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FIG. 3: The predicted first-order p3| octupole deformations 
(top), intrinsic dipoie moments (middie) and intrinsic SchifF 
moments (bottom) for four Skyrme interactions in a series of 
even-A^ Radium isotopes. The experimental intrinsic dipoie 
moments are also shown. Where symbols are missing, the cor- 
responding predicted values are zero because the mean field 
is not asymmetrically deformed. 



given in Ref. 30] . (In ^^"'Ra and ^^^Ra, by the way, 
we agree fairly well with the quadrupole moments ob- 
tained from B{E2ys in Ref. For example, SkO' gives 
(32 = 0.184 in ^^'^Ra and experiment gives /§2 = 0.179(4).) 

The second panel in the figure shows the absolute val- 
ues of intrinsic dipoie moments Dq = e (zp) , along 
with experimental data extracted from El transition 
probabilities [22|. The calculated values for Dq change 
sign from positive to negative between N = 134 and 
N ~ 138, reflecting a small change in the location of the 
center of charge from the "top" half of the pear-shaped 
nucleus to the "bottom" half. This predicted sign change 
is consistent with the near-zero experimental value for 
N = 136. None of the forces precisely reproduces the 
trend through all the isotopes, but the comparison has 
to be taken with a grain of salt because "data" derive 
from transitions between excited rotational states, and 
therefore are not necessarily identical to the ground-state 
dipoie moments. Cranked Skyrme-HF calculations with- 
out pairing correlations [2J|and cranked HFB calcula- 
tions with the Gogny force [23 predict that for most Ra 
isotopes Dq changes significantly with angular momen- 
tum. In any event, as thoroughly discussed in Ref. [2^ . 
the intrinsic dipoie moment is a small and delicate quan- 
tity. 

The intrinsic SchifF moment (Sz), the quantity we're re- 



ally interested in, is more collective and under better con- 
trol, as the bottom panel of the figure shows. The various 
predictions are usually within 20% of one another and 
large, confirming the predictions originally made in Rcfs. 
[^Q. The octupole deformation and intrinsic dipoie mo- 
ment have been shown to change onl y sl ightly with parity 
projection from the intrinsic states ^3], and the same is 
probably true of the intrinsic Schiff moment. 

By turning the pairing force off, we are able to 
see whether the parity-violating quantities in Fig. |3| 
are affected by pairing correlations. In ^^^Ra, for 
example, SkO' gives /Jg = 0.141, Dq = -0.103 efm, 
and {Sz) = 34.4 efm"^ without pairing, and (3^ — 0.143, 
Dq = —0.093 efm, and {Sz) = 34.3 efm^ when pairing is 
included. In this nucleus uncertainties related to pairing 
are very small. 




FIG. 4: Single-particle spectra for protons (top) and neutrons 
(bottom) in ^^''Ra, for the four Skyrme interactions. 

Finally, in Fig. 21 we show the predicted proton 
and neutron single-particle spectra generated by the 
ground-state mean-field in ^^"'Ra. The combination of 
quadrupole, octupole, and higher deformations reduces 
the level density around the Fermi surface for both kinds 
of nucleon, leading to significant deformed Z — 88 and 
N — 132 shell closures for all interactions, and a some- 
what weaker N = 136 subshell closure for SIII, SkM* 
and SkO'. The small level density around the Fermi sur- 
face might explain the insensitivity of the deformation to 
pairing correlations mentioned above. For all the forces 
except SkM*, the first empty neutron level clearly has 
jz — 1/2, implying that in ^^^Ra the ground-state parity- 
doublet bands will be built on J'^ = 1/2* states. For 
SkM* the situation is less clear because the jz = 1/2 and 
3/2 states are nearly degenerate, and it is necessary to 
carry out the calculation in ^^^Ra itself to see which be- 
comes the lowest. 
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IV. CALCULATING THE SCHIFF MOMENT 

OF ^^"Ra 

A. Odd-yl Nuclei and SchifF Moments in 
Mean-Field Approximation 

Fully self-consistent calculations in odd-^ nuclei are 
possible but seldom performed. For many physical ob- 
servables it is enough to neglect correlations between the 
odd particle and the core, which amounts to dropping 
a valence particle into the field generated by that core. 
Other quantities, however, are sensitive to the interac- 
tion between the last particle and the core. The inter- 
action can change the deformation and pairing strength, 
and produce various kinds of core polarization. In fact, 
a self-consistent odd-A calculation is equivalent to first 
performing the calculation in the even-even nucleus with 
one less neutron, then placing the last nucleon in the first 
empty orbit and treating its polarizing effect on the core 
in RPA [33, HI]. Fully self-consistent calculations take all 
kinds of polarization into account simultaneously. For us, 
the spin polarization is of particular interest. 

To proceed here we use the public-domain Skyrme- 
Hartree-Fock code HFODD (v2.04b) [Ulllll^, which 
allows the mean-field and the associated Slater determi- 
nant to simultaneously break invariance under rotation, 
reflection through three planes, and time reversal (the 
code we used in the last section does not). Breaking 
the first is necessary to describe a deformed nucleus. 
Breaking all three reflections (and not only parity) is 
necessary to represent axial octupole deformation with 
the spin aligned along the symmetry axis. Breaking the 
last induces spin polarization in the core, which because 
of Kramers degeneracy cannot otherwise occur. Incor- 
porating spin polarization is important because it has 
the potential to signiflcantly alter the matrix element of 
tjpT in Eq. (jSJ from its value in the particle-rotor model, 
where the spin is carried entirely by one valence parti- 
cle. The code HFODD cannot yet treat pairing when it 
allows T to be broken, but pairing in T-odd channels is 
poorly understood. No existing codes can do more than 
HFODD in odd-A octupole-deformed nuclei. 

As above, we use the Skyrme interactions SIII, SkM*, 
SLy4, and SkO'. The reason SkO' is our favorite has to 
do with the part of the energy functional composed of T- 
odd spin densities (which, following common practice, we 
refer to as the "T-odd functional" , even though the entire 



functional must be even under T) . The T-odd functional 
plays no role in the mean-fleld ground states of even nu- 
clei, but can be important in any state with nonzero an- 
gular momentum. Of the forces above, only SkO' has 
been seriously investigated in T-odd channels. In Ref. 
[TtII . the T-odd part of the functional was adjusted to re- 
produce Gamow- Teller resonances, resulting in an effec- 
tive Landau parameter (/g = 1.2. In the isoscalar channel, 
the force was adjusted to reproduce the commonly used 
value go = 0.4 [33]. Although there are not enough data 
to constrain other relevant parameters in the functional, 
and although a very recent calculation starting from a 
realistic nucleon- nucleon interaction 38], while confirm- 
ing the value g'^ — 1.2, finds go = 0.85, SkO' is clearly 
the best available Skyrme interaction for describing spin- 
spin interactions. The corresponding T-odd terms in the 
functional are precisely those that will polarize the spin 
in the core. There are other terms that could be added to 
the standard Skyrme interaction and do the same thing 
— - the tensor force for example — but they are almost 
never used and their effects still need to be investigated. 

The parameters of the other three forces SIII, SkM*, 
and SLy4 were adjusted entirely to ground-state proper- 
ties in even-even nuclei, and so the Landau parameters go 
and go were not fit. Here we set them to zero by treating 
the T-odd and T-even terms in the Skyrme functional 
independently, as described in Ref. [T^l- Only orbital 
terms, which are fixed by gauge invariance |33 (a gen- 
eralization of Galilean invariance), appear in the T-odd 
parts of these forces. 

We rely on the crudest forms of projection. For parity, 
that means none at all, and for angular momentum it 
means inserting the rigid-rotor factor J/{J + 1) = 1/3 in 
front of the intrinsic Schiff moment, as described above. 
In other words, we use Eq. lO, with the intrinsic state 
taken to be the Hartree-Fock ground state produced by 
HFODD. Just as in the particle-rotor model, the intrinsic 
Schiff moment is given by the classical expression Eq. J^l, 
but with pch equal to e times the Hartree-Fock ground- 
state proton density. As already mentioned, the Hartree- 
Fock approximation allows us, by summing over occupied 
orbits, to write the intrinsic matrix element of the two- 
body potential VpT as the expectation value of an ef- 
fective one-body operator UpT- Because we now have a 
microscopic version of the "core" , this effective potential 
is more complicated than in Eq. (jS)), and it now acts on 
all the nucleons: 



PT 



A 

1=1 



• V 



d^ 



47r|r - 



(.90 + 2g2)pi(r') - gipoir') 



exch. 



(9) 



Here po{r) = Pn{r)+ Pp{r) and pi{r) = pn[r)~ pp[r) are the isoscalar and isovector densities. The piece coming 
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from exchange terms in the original two-body interaction 
VpT is nonlocal, just as in the usual Hartree-Fock mean 
field, and we have not written it explicitly here (though 
we do below). 

The code HFODD at present cannot evaluate the ex- 
pectation value of a folded potential like that above, 
which is due to the finite pion range. Nevertheless, even 
in the zero-range approximation we can avoid several of 



the assumptions — proportionality of neutron and proton 
densities, negligibility of exchange terms, and absence of 
core spin — leading to the extremely simplified potential 
in Eq. (jSJ . The zero-range approximation is equivalent to 
assuming the pion is very heavy, so that the term involv- 
ing the pion mass in Eq. (j^ becomes a delta function. 
Under that assumption, but none others, the exchange 
terms become local and UpT takes the form: 



UpTir) 



1=1 
1 ^ 

+ 2^^"-' i^^9o + giTz,t)Jo{r) + {gi + goTz.i - 'i92Tz,i)Ji{r) |. 



(10) 



Here J( r) is the "spin-orbit" current, defined, e.g., in 
Ref. and references therein, and the subscripts and 
1 refer to isoscalar and isovector combinations as they 
do for the density. The terms in UpT that contain J 
are the exchange terms omitted above. We will evaluate 
them, but argue later that their effects are probably small 
when the finite range is restored. The terms containing 
the density p all come from the direct part of VpT- We 
do not simplify things further to obtain something like 
Eq. © because pp is not really proportional to p„ and 
the core nucleons do carry some spin. We will manage 
nevertheless, to compare our results with those of Ref. 
@ . We will also estimate the effect of a finite pion range 
on the direct terms, though our inability to do so more 
precisely at present is the most significant shortcoming 
of this work. 

HFODD works by diagonalizing the interaction in the 
eigenbasis of an optimal anisotropic three-dimensional 
harmonic oscillator. For ^^^Ra, algorithms developed in 
Ref. 34] give oscillator frequencies of hujz— 7.0625 and 
huj±=8.6765 MeV in the directions parallel and perpen- 
dicular to the elongation axis. The matrix element of 
UpT converges only slowly as we increase the number of 
levels in the basis. When the interaction polarizes the 
core, it takes 2500 or more single-particle basis states to 
get convergence. The basis then contains up to Nz— 26 
and N±=21 oscillator quanta. 

B. Laboratory Schiff Moment of ^^'^Ra 

We turn finally to results in ^^^Ra itself. For SkO', our 
HFODD calculations yield P2 = 0.190, P3 = 0.146, and 
(34 = 0.136 for the usual first order approximation to the 
deformation parameters determined from mass multipole 
moments 29]. The laboratory Schiff moment, Eq. 10, is 
proportional to the product of the intrinsic Schiff moment 
(Sz) and the expectation value (Vpt)- Table ^ shows 



the intrinsic moments and the expectation values of the 
6 operators that enter the zero-range approximation to 
VpT in Eq. H10|l . Before commenting on the entries, we 
mention what is in each of the forces and calculations. 

For all the forces, terms in the functional that are pro- 
portional to Laplacians of spin densities (s • As) and 
density-dependent spin-spin terms {f{p)s ■ s), cf. Ref. 
jirl [33, which enter through the T-odd part of the 
Skyrme functional, have been turned off. For the first 
three lines in Table HI [forces labeled with (0)] , the spin- 
spin terms have also been turned off, so that the only 
nonzero terms in the T-odd functional (as noted above) 
are those required by gauge invariance. For the fourth 
line [SkO'( — )], all T-odd terms in the functional have 
been turned off, so that aside from the self-consistency 
in the wave functions the calculation resembles one with 
a phenomenological (non-self-consistent) potential, for 
which T-odd mean-fields are never considered. We in- 
clude this result so that we can distinguish the role 
played by core polarization. The results in the line la- 
beled SkO' include the time-odd channels, adjusted as 
discussed above ^3- This is the force in which we have 
the most confidence. The last entry is the result of Ref. 

with the implicit assumption that the neutron and 
proton densities are proportional. 

In our calculations, the intrinsic Schiff moments are 
close to one another, and all are less than twice the es- 
timate of Ref. . The agreement reflects the collective 
nature of these intrinsic moments; they are even larger 
than the particle-rotor estimates. But the matrix ele- 
ments of VpT, the other ingredient in Eq. ^ for the 
laboratory Schiff moment, are a bit more delicate. Our 
results show the exchange terms on the right side of the 
table to be comparable to the direct terms, a result that 
is surprising because for a spin-saturated core (or in the 
particle-rotor model) the exchange terms vanish exactly. 
We think, however, that the ratio of exchange to direct 
terms would become small were the finite range of the in- 
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TABLE I: The intrinsic SchifF moment, in units of efm'^ and the intrinisic-state expectation values of operators in Eq. IIUH . in 
units of lO'^fm"*. 





(o-r ■ Vpo) 


(o-r ■ Vpi) 


{<T ■ Jo) 


(a . Ji) 


(o-r • Jo) 


(o-r- Ji) 


SIII(O) 


34.6 


-1.081 


-0.087 


-1.047 


0.162 


-1.049 


0.159 


SkM*(0) 


46.6 


-0.730 


-0.497 


-1.043 


0.099 


-1.042 


0.085 


SLy4(0) 


43.4 


-0.676 


-0.578 


-1.303 


-0.016 


-1.299 


-0.019 


SkO'(— ) 


41.7 


-0.343 


-0.318 


-1.149 


0.030 


-1.149 


0.030 


SkO' 


41.7 


-0.467 


-0.227 


-0.476 


0.014 


-0.526 


0.014 


Rof. [6] 


24 


-2 


-0.4 











teraction reintroduced and short-range NN correlations 
inserted. 

Though unable to include either effect here, we did so 
in a Nilsson model for ^^^Ra. We took nucleons there 
to occupy independent single-particle levels generated by 
a deformed oscillator potential with (32 = 0.138, (3^ = 
0.104, and /34 = 0.078, values taken from Ref. ^J. We 
then evaluated the ground-state expectation value of the 
full two-body interaction VpT, with and without the zero- 
range approximation (and in the latter case, with short- 
range correlations included a la Ref. 40] ). In this simple 
model, the valence nucleon carries all the spin, and only 
the neutron-proton and neutron-neutron parts of VpT 
contribute. The direct np term shrank by a factor of 
1.5, while the corresponding exchange term shrank by a 
factor of 1400 (both independently of the ^'s in Eq. Q), 
it turns out) when the range of the interaction was set 
to its proper value. The results in the nn channel were 
less dramatic: the direct part again shrank by 1.5 and 
the exchange part by a factor of 5. When we moved the 
valence neutron to higher orbits, these numbers changed 
some — the direct terms sometimes were not suppressed 
at all and other times shrank by factors of up to 6, but 
the ratios of the exchange to direct contributions almost 
always ended up small. Similar behavior was found for 
parity- violating forces in Ref. , where it was traced in 
part to the different average momenta carried by the pion 
in direct and exchange graphs. So that we can compare 
our results with those of Ref. pf, we will neglect the 
exchange terms from now on, though we caution that 
this step should eventually be justified more rigorously, 
e.g., by actually calculating them with the finite-range 
force in the full mean- field model. The reduction we see 
in the direct terms is in line with the results of Ref. j 
though we find it more variable^. 

Though we cannot yet be more quantitative about 
finite-range effects, we do quantify the core polarization 
in Table For the first three lines of the table, where 



^ We performed another test, using the direct part of Eq. IIUI 
with the valence wave function taken from the Nilsson model just 
described, but with the neutron and protons densities assumed 
to have more realistic Wood-Saxon forms. The direct terms were 
again suppressed by factors of 1.5 to almost 10 that depended 
significantly on the valence orbit. 



the forces are labeled (0), the spin-spin terms are ab- 
sent from the energy functional, and the protons in the 
core develop only a tiny spin density from the T-odd 
terms required by gauge invariance. For the fourth line, 
SkO'( — ), all T-odd terms are absent and the protons 
can have no spin at all. This means that the operators 
f{r)cr and f{r)(TT have either the same or almost the 
same expectation value for any /(r) so that columns 4 
and 6 ((<t- Jq) and {(tt-Jq)) have identical or nearly iden- 
tical entries for these forces, and so do columns 5 and 7 
(((T • Ji) and (err • Ji)). The fifth line of the table con- 
tains the effects of spin polarization, which are primarily 
to alter the neutron-spin density; the equalities between 
the columns are not badly broken, so the protons do not 
develop much spin. The same is true of the terms involv- 
ing p, though that is not obvious from the table because 
we display only the two terms that appear in Eq. pO(l . 

These near equalities and the probable irrelevance of 
the exchange terms when the finite range is taken into 
account imply that only the quantities cr„ • V/9„ and 
cr„ • Vpp are ultimately important. We display them 
in Table ITU Except for SIII, the neutron-density distri- 
bution affects the matrix element much more than the 
that of protons. By comparing the fourth and fifth lines, 
however, we see that spin correlations increase the role of 
the protons, while reducing that of the neutrons slightly. 
Thus, while the spin-spin interactions do not cause the 
protons to develop much net spin, they do correlate the 
neutron spin with the proton density. 



TABLE II: Intrinsic-state expectation values of important 
matrix elements, in the neutron-proton scheme, in units of 
10-^fm-*. 





(o-„ • Vp„) 


(o-„ ■ Vpp) 


SIII(O) 


-0.577 


-0.491 


SkM*(0) 


-0.619 


-0.120 


SLy4(0) 


-0.628 


-0.050 


SkO'(— ) 


-0.331 


-0.013 


SkO' 


-0.320 


-0.114 


Ref. f6] 


-1.2 


-0.8 



There is not too much scatter in our results. The en- 
tries in the second column ((<Tr • Vpo)) of Tablcnjdiffer 
by factors of two or three, and the entries in the third 
((err • Vpi)) by a little more, though they are all smaller 
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than those in the second column (which is not surprising 
— the third column subtracts the neutron and proton 
densities while the second adds them). In the neutron- 
proton scheme (tabled all of our numbers are smaller 
than those of Ref. ^J, a result that was anticipated in 
Ref. 0. The difference from the earlier estimate for the 
larger matrix elements ranges from factors of two to four, 
though the isovector combination — the third column in 
table U — is sometimes actually enhanced a little. 

What, at last, have we to say about the real labora- 
tory Schiff moment 5"? The lab moment is given by the 
product of the matrix elements just discussed, the in- 
trinsic Schiff moments, and the unknown coefficients gi. 
Our intrinsic Schiff moments are about 1.5 times larger 
than those of Ref. while our VpT matrix elements, 
in the zero-range approximation, are smaller than theirs, 
usually by a somewhat larger amount. Overall, our lab 
moments will usually be smaller by factors between about 
1.5 and 3 than the estimates of Ref. 6] (an exception can 
occur if for some reason gi is considerably less than the 
other two coefficients). 

How big are our moments compared to that of ^^^Hg? 
The most comprehensive calculation in that nucleus, 
which appeared very recently 43j , improved on the work 
of Ref. Q by including the effects of the residual strong 
interaction and the full finite-range form for VpT- The 
new results are smaller than that of ref. , only slightly 
so for the isovector part of VpT, but by a considerably 
amount in the isoscalar and isotensor channels. The au- 
thors write their results in terms of the pion-nucleon cou- 
plings as 

S'Hg = .0004 ggo + .055 ggi + .009 gg2 [efm^]. (11) 

Our result for radium, with the zero-range approximation 
and exchange terms neglected, translates to 

^^ero-rangc ^ _5 ^QQ ggo+W A gg.-lO .1 gg, [efm\ (12) 

If the three g's are comparable, therefore, our Schiff mo- 
ment is several hundred times larger than that of Ref. 
[43l|, in part because the isoscalar and isotensor interac- 
tions are more effective in Ra than in Hg. [If gi is larger 
than the other two couplings, as in left-right symmet- 
ric models [i^, our result is less than 200 times bigger 
than the latest one in ^^^Hg. The very small coefficient 
of go for ^^^Hg in Eq. Ijllll . by the way, has significant 
consequences |43 for the limit on the QCD T-violating 
parameter that can be inferred from the experimen- 
tal limit in Ref. ^ .] Accepting the work of Ref. on 
atomic physics in Ra and Hg, the enhancement of the 
atomic EDM of ^^^Ra is about about three times that 
of the Schiff moment, i.e. potentially more than 1000. 
We again caution, however, that we have yet to include 
the full finite-range version of Vpx and that our prelim- 
inary investigations suggest that doing so will decrease 
our Schiff moment at least a little. Ironically, Ref. 
finds that including the finite range actually increases the 
matrix element in ^^^Hg, though only slightly. 



We hope to make other improvements in our calcu- 
lation as well. Projection onto states of good parity 
will change the results a bit, as will angular-momentum 
projection. Our conclusions about the size of spin- 
polarization corrections could be modified by two terms 
in the Skyrme functional we've set to zero, or by a better 
determined value of the Landau parameter go. Whatever 
the result of such corrections, however, it is clear that the 
atomic EDM of ^^^Ra will always be significanlty larger 
than that of ^^^Hg. 



V. CONCLUSIONS 

We have calulated the Schiff moment in ^^^Ra in the 
approximation that the T-violating interaction has zero 
range. Our calculations, which are self-consistent and 
include core polarization, give results that are generally 
just a few times smaller than earlier estimates based on 
the particle-rotor model. Accepting the very recent re- 
sults of Ref. ^3 , we currently find the Schiff moment of 
^^^Ra to be (generically) several hundred times that of 
i99jjg^ a result that strengthens the case for an atomic- 
EDM experiment in Ra, though the enhancement factor 
depends significantly on the source of T violation, and 
we expect it to decrease at least a little when we use 
the finite-range force. Work towards including a finite 
range in HFODD is in progress. We also plan to apply 
the self-consistent methods used here to other light ac- 
tinides, as well as to ^^^Hg, where we suspect octupole 
correlations may play some role Q . Maintaining self con- 
sistency in i^^Hg should automatically control the spuri- 
ous Schiff strength encountered in Ref. . The source 
of the insensitivity of the Schiff moment to T violation 
in the isoscalar channel in that work should be checked 
and understood. 

After many years of neglect, the question of which 
isotopes are best for EDM measurements is now being 
rapidly addressed. 
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